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The concept of collective [lj degrees of freedom has 
been of central importance in nuclear physics, if only to 
generate models involving much fewer degrees of free- 
dom than the true number, 3A, needed for a proper mi- 
croscopic description of a nucleus of mass number A. 
Under adiabatic assumptions, or even proofs, the de- 
scription of nuclear dynamics can be compressed into 
slow motions of a few collective degrees of freedom B, 
while the other degrees, faster, can be averaged out. Of- 
ten, such collective degrees can be one-body operators, 
B = ^ i=1 /?(ri,pi,a-i,Ti), where r^p^a^n refer to the 
position, momentum, spin, and isospin, respectively, of 
nucleon i. The summation over all nucleons provides, in- 
tuitively at least, a motivation for more inertia in B than 
in the individual (nucleonic) degrees of freedom /3j. 

The concept of energy surfaces Q has been almost as 
important. Given a "coordinate- like" collective opera- 
tor B and its expectation value b = (B), every collec- 
tive model accepts the common wisdom that there exists 
an energy function, e(b), and an inertia parameter, /x(6), 
that drive the collective dynamics. This concept is at the 
core of any theory of fission or fusion, for instance. Key- 
words like "saddles", "valleys", "barriers" , etc., flourish. 

Simultaneously, it is implicitly assumed that the func- 
tion, e(6), results from an energy minimization under 
constraint. Namely, while the system evolves through 
various values of b, it is believed to tune its energy to 
achieve a local minimum. To illustrate, consider the nu- 
clear Hamiltonian, H = ^ ti + J2i<j v ij> where t and v 
denote the usual kinetic and interaction operators. Given 
a suitable trial set of density operators T> in many-body 
space, normalized by TtD = 1, a theory for an energy 
surface most often reads, 



e(b) = inf Tr {HV\ 

V^b 



(1) 



where Tr is a trace in the many-body space for A nucle- 
ons and the constraint, V b, enforces Tr {BT>} = b. 
In practice, this requires a diagonalization, or at least a 
minimization, of the constrained operator, % = H — XB, 
where A is a Lagrange multiplier. The purpose of this 



work is to show that this constrained minimization can 
only return concave 0] functions e(b). Maxima are 
impossible. In the trivial generalization where several 
collective operators B\, . . . , B/v, are involved, concavity 
stills holds, hence saddle points are excluded as well as 
maxima, whether local or not. Only an absolute mini- 
mum is available. 

Consider, indeed, the exact diagonalization at T = 0, 
or an exact partition function at finite T. However, if A is 
a positive number and the spectrum of B is not bounded 
from above, there is a risk that the spectrum of T~L extends 
to — oo, hence no ground state of H exists to provide a 
minimal energy. A similar risk occurs if A is negative and 
B is not bounded from below. To start, we shall assume 
that either H is always bounded from below or both H 
and B have been projected in a finite subspace of A- 
body basis states; in the latter case, the problem reduces 
to handling finite matrices. Let ip(\) be the ground state 
of H. (For the sake of simplicity, we assume that there 
is no degeneracy.) The corresponding eigenvalue, e(A), is 
stationary with respect to variations of ip, among which is 
the "online" variation, dX (dip/dX), leading to the well- 
known first derivative, de/dX = —b = — (ip \B\tp) . Let 
P = \ip) (i/>\ and Q = 1 — P be the projector and the 
complementary projector, respectively. Brillouin-Wigner 
perturbation theory yields the first derivative of ip, viz. 
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Since the operator, (e — QHQ), is clearly negative- 
definite, the energy, e, is a convex function of A. It 
is trivial to prove that the same convexity holds for 
the ground state eigenvalue e(Ai, . . . , Xn) if several con- 
straints, B\, . . . , B/v, are used. When a temperature T 
is considered, a partition function is calculated from an 
H bounded from below, and one obtains a free energy, 
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e(Ai, . . . , Ajv;T), that also contains the entropy contri- 
bution, -TS, where S = -Tr{V\nV}. A proof of the 
convexity of an exact e(Ai, . . . , X^;T) is also easy 

At T = 0, the traditional Legendre transform ex- 
presses the energy, e = (ijj \H \ ip), in terms of the con- 
straint value(s) rather than the Lagrange multiplier (s). 
For simplicity, set one constraint only, the generalization 
to N > 1 being easy. Since e = e + Xb, and de = — bdX, 
then de = Xdb, giving de/db = X, a most traditional re- 
sult for conjugate variables. Furthermore, the second 
derivative, d 2 e/db 2 , reads, dX/db — \/{db/dX). From 
Eq. ([3]), the derivative, db/dX, is positive-definite. Ac- 
cordingly, e is a concave function of b. Now, if T > 0, 
the Legendre transform instead generates a reduced free 
energy, r/ = (e — TS 1 ) , which is a concave function of the 
constraint value(s). An additional Legendre transform, 
between T and S, returns e alone, as a concave function 
of both the constraint(s) and S. 

Let b- and 6+ be the lowest and highest eigenvalues of 
B, respectively. It is expected that b spans the interval, 
[&_,&+], when A runs from — oo to +oo. In every exact 
diagonalization of W, or exact partition function calcula- 
tion, concavity sets a one-to-one mapping between b in 
this interval and A. More generally, with exact calcula- 
tions, there is a one-to-one mapping between the set of 
Lagrange multipliers, {Ai, . . . , Ajv}, and that of obtained 
values, . . . , 6jv}, of the constraints. Concavity, in 
the whole obtained domain of constraint values, imposes 
a very poor landscape: one valley only. 

Consider, now, situations where the energy surface 
does not result from a diagonalization, or the exact cal- 
culation of a partition function. Typically, the min- 
imization could be estimated from a Hartree-Fock(- 
Bogoliubov) approximation, whether at zero or finite T. 
The trial states involved in such nonlinear approxima- 
tions do not correspond to a linear manifold; indeed, for 
instance, a sum of two determinants usually does not 
make a determinant. Let T>(X) denote one A-body den- 
sity operator where, within such nonlinear approxima- 
tions, a minimum of Tr {HT>} or, at finite T, a minimum 
of (Tr {HV} — TS) is reached. Let e(A) denote this min- 
imum value. It may be degenerate, incidentally, because 
of symmetry breaking but, whether degenerate or not, 
this value is stationary with respect to arbitrary varia- 
tions 5T> within the manifold of trial states. As a con- 
sequence of this stationarity, the first derivative again 
reads, de/dX = — b = — Tr {BT>}. Then, if a Legendre 
transform is possible, defining r\ = e + Xb, in terms of 
b, the same argument that was used for the exact case 
again yields, dr\jdb = X. With several constraints, the 
gradient of 77 in the vector space spanned by {61, . . . , bjsr} 
is again the vector {Ai, . . . , Xn}- 

There remains to discuss second derivatives. Con- 
sider, for instance, Hartree-Fock (HF) calculations, where 
A-body density operators are dyadics of Slater deter- 
minants, T> — \tp) (ip\. Norm-conserving variations of 



an HF solution, <p, can be parametrized as, \8<p) — 
exp(iXSa) \<p) , with A an arbitrary particle-hole Her- 
mitian operator, and 5a an infinitesimal coefficient. Un- 
der such a variation in the neighborhood of an HF solu- 
tion, the first and second order variations of the energy, 
e = Tr {Hexp(iXSa)T>exp(—iXSa)}, read, 

Se = i5aTr{[n, X]V} = 0, (4) 

and 

S 2 e = - (8a 2 /2) Tr{[[% X] ,X]V}> 0, (5) 

respectively. Since we are dealing with an HF solution, 
the first order vanishes VA, and since only those solutions 
that give minima are retained, the second order variation 
of e is semi-positive-definite, VA again. Now, when % 
receives the variation, —BdX, there exists a particle-hole 
operator Y, a special value of A, that, with a coeffi- 
cient dX, accounts for the modification of the HF solu- 
tion. This reads |$) = exp(iy<iA) \<p). Simultaneously, 
those particle-hole operators that refer to this new Slater 
determinant $ become X = cxp(iYdX)X exp(-iYdX). 
The new energy is thus, 

£ = Tr {exp(-iYdX)(n - BdX) cxp(iYdX)V} . (6) 

The stationary condition, Eq. (j4}, becomes, 

= Tr {exp{-iYdX) [(H - BdX),X] exp{iYdX)V} 
= Tr {[exp(-iYdX)(H - BdX) cxp(iYdX), A] V} . (7) 

The zeroth order in dX of this, Eq. 0, reads, 
Tr{[H,X]V}, and vanishes VA, because of Eq. ©. 
Then the first order in dX gives, again VA, 

Ti{[B,X]V} = iTx{[[n,Y],X]V}. (8) 

The second derivative is, 

d 2 e/dX 2 = -(d/dX)Ti{exp(~iYdX)Bexp(iYdX)V} 

= -tTr{[B,Y]V}. (9) 

Upon taking advantage of Eq. (jHJ, for Y as a special case 
of A, this becomes, 

d 2 e/dX 2 =Tr{[[H,Y],Y]T>}, (10) 

the right-hand side of which is semi-negative-definite, see 
Eq. (|5|). The solution branch obtained when A runs is, 
therefore, convex. Its Legendre transform is concave. 

Beyond HF, the proof can be generalized as follows. 
Consider a solution branch X'(A), and expand D(X) up 
to second order, assuming that the manifold of solutions 
is suitably analytic, 

V(X + dX) = D(A) + dX{dV/dX) + (dX 2 /2)(d 2 V/dX 2 ). 

(11) 
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Let A receive a variation dX in V(X) alone, not in H. The 
stationarity and minimality of Tr {HT>} with respect to 
any variation of V induce, 

Tr {ndV/dX} = 0, 
Tr {Hd 2 V/dX 2 } > 0. (12) 

There is also a stationary property for the solution, 
V(X + dX), but the Hamiltonian is now, H(X) — BdX, 
and the derivative of the state becomes, dV/dX + 
dX(d 2 V/dX 2 ) + 0{dX 2 ), and therefore, 

= Tr {(% - BdX) [dV/dX + dX(d 2 V/dX 2 ) + 0{dX 2 )] } . 

(13) 

The zeroth order of this, Eq. (13]), reads, Tr {HdV/dX}. 
It vanishes, because of the first of Eqs. (fT2|) . The first 
order, once divided by dX, gives, 

- Tr {BdV/dX} = -Tr{nd 2 V/dX 2 } . (14) 

The left-hand side of Eq. (jT4"]) is nothing but the the 
second derivative, d 2 e/dX 2 . The right-hand side is semi- 
negative-definite, because of the second of Eqs. (TT21 . 
Concavity of e(b) is again a property of such approxi- 
mations with strict minimizations. 

A landscape can be recovered, however, in a some- 
what contrived way, if one accepts the crossing of var- 
ious solution branches. Exact diagonalization normally 
prevents level crossing, unless distinct quantum numbers 
are considered, but approximations, usually nonlinear, 
often generate multiple solutions corresponding to local 
minima. Special solutions may also occur if, in special 
cases, variational spaces have edges because the energy 
gradient does not need to vanish: it just has to be or- 
thogonal to the edge. We assume, without proof, that 
such exceptional edge solutions occur either only for iso- 
lated values of A or, if they make true functions e(A), that 
they retain enough analyticity for a concavity of the cor- 
responding e(b). All told, various branches of e(b) may 
intersect, and the pattern formed by their lowest parts 
is, piecewise, concave. This is schematically illustrated 
by Figs. [T]and [2] Fig. [T] shows the plots of three convex 
parabolas, with equations, e = — 2 — 6A — A 2 , e = — 1 — A 2 , 
and e = —2 + 4A — A 2 , respectively. Then Fig. [5] shows, 
from their Legendre transforms, the plots of e(b), namely 
e = 7 + 36 + 6 2 /4, e = -1 + & 2 /4, e = 2 - 2b + 6 2 /4, re- 
spectively. Their lowest segments make a pattern with 
interesting structures. The physical significance of the 
spikes, however, is doubtful. A further theory for the 
rounding of such singularities seems to be in demand. 
A calculation of collective inertia parameters smoothly 
through them does not seem to be obvious either. One 
would have to face the contradiction with the "one valley 
syndrome" of an exact solution. How could an approxi- 
mation yield better physics than the exact solution? 

One may argue that the zoo of stationary solutions 
of such approximate methods is usually rich enough to 




FIG. 1: Toy example with three branches e(A), as discussed 
in the text. 
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FIG. 2: Legendre transforms, e(6), of the branches in Fig. Q] 

accommodate inflections, where the curvature signature, 
Tr {'Hd 2 'D/dX 2 } , can change sign. "Phase transitions", 
a somewhat incorrect wording for a finite system, are 
sometimes advocated to accept continuing branches of 
energy minima into metastable branches. But this means 
dropping the hypothesis of strict energy minimization. 

For unbounded constraint operators, a further paradox 
may occur, as seen from the following example. Assume 
that the interaction v present in H has a strictly finite 
range, namely Vij = if |rj — ij\ > r$. Consider a set 
of orbitals, tpi, . . . , ipA, with wave functions (fi(r, a, r) = 
Cj exp [— l/(r - Ri) - 1/(S< - r)] /r if R% < r < S i: and 
(fi = otherwise. Here r = |r|, and Cj is the normal- 
ization coefficient. Given a large positive number L, 
set Si = R\ + L, . . . , Sa = Ra + L, so that such or- 
bitals are very flat, and, therefore, carry a minuscule 
kinetic energy, cx 1/L 2 , as small as wanted. Set, also, 
R2 > Si + ro, ... , Ra > Sa-i +fo, so that such orbitals 
are so separated that, not only their orthogonality is triv- 
ial, but also any antisymmetrized interaction matrix ele- 
ment, {tpiPj \v\ (pk<pe), vanishes. The Slater determinant 
X made of such orbitals has an energy, (x\H\x), as close 
to zero as wanted. Simultaneously, for a collective oper- 
ator such as, B v — Ylt=i r ii with, typically, v > 2, it is 
trivial to show that AR\ < (x {B^x) < AS A . The situa- 
tion is illustrated in Fig.|3]for a chosen set of parameters. 
Now, given the ground state of H, with its eigenvalue 



4 



0.05 




0.04 




0.03 




0.02 




0.01 









20 40 60 80 100 



FIG. 3: Example, with i?i = 20 fm, L = 23 fm and r = 2 fm, 
see text, of orbitals ifi for a Slater determinant \ that, if 
mixed with the ground state, can leave the energy hardly 
changed while allowing a very large collective moment. True 
ground state nuclear densities, halos and tails included, oc- 
cupy the first 10 fm at the very most. 

E , define an A-body state by the normalized mixture of 
density matrices, V mix = (|* > (*o X) (xl) /(I 
where the positive number u> will eventually become in- 
finitesimal. When L — > +oo, the energy of this mix- 
ture can be made as close to E /(l + lj) as wanted. 
Simultaneously, if R\ — > +oo, the expectation value, 
b v = Tr {B v V m i X }, of the collective operator diverges 
at least as fast as ujAR"/{1 + uS). Maintain constant the 
product, loRi, and the ratio, L/Ri, while R\ grows to in- 
finity. It is clear that the limit of the energy is Eg , while 
b v can grow arbitrarily, like oc R\~ x at least. This makes 
the right-hand side of the energy surface, e(b v ), flat. 

The paradox still holds with interactions v that do not 
strictly vanish beyond tq. For any tail of v, the separa- 
tions, R% — Si, . . ■ , Ra — Sa-i can be adjusted to make 
the matrix elements, {fi^Pj \v\ ipkPi), as small as wanted. 

For traditional nuclear multipole operators, a spherical 
harmonic multiplies the radial factor, r v . We insert the 
same harmonic in the </?'s to extend the paradox; an arbi- 
trary value of a multipole is compatible with the ground 
state energy. A choice of collective operators that are 
bounded seems to be mandatory if the concept of a non- 
flat landscape for a collective energy has to be salvaged. 

Even under such a precaution, this work shows that 
concavity, piecewise at least, is a major property of any 
energy surface obtained by a strict minimization of the 
energy under constraint (s). If an energy landscape with 
"mountains" and "saddles" is needed, this concavity con- 
tradicts the intuition that energy transits through local 
minima. The success of collective models that use a non 
trivial landscape is too strong to be rejected as physically 
and/or mathematically unsound, but its validation must 
now be searched through other methods such as, likely, 



resonating group methods generator coordinate ones 
Born-Oppenheimer approximations, influence func- 
tionals Q, deconvolutions of wave packets in collective 
coordinate spaces 0], etc. Except for anharmonic vibra- 
tions, where one valley is sufficient, one might have to 
accept that metastability may be as important as en- 
ergy minimization. Or could it happen that collective 
operators have to be A-dependent, in order to convert 
metastable branches for constant operators into strict 
minimization branches for A-dependent ones? Since fluc- 
tuations are important at "phase transitions" , combina- 
tions of the form, K, = H — AB — f(A)B 2 , with at first very 
simple forms of /(A), deserve some attention. Such oper- 
ators /C govern both a constraint and its fluctuation, but 
obviously differ from a double constraint form, H = H — 
X\B — A2-B 2 , with two independent parameters, Ai,A2. 
A one-dimensional path with non trivial structures, be- 
cause of an additional term, -d 2 f/d\ 2 ({B + B 2 df/dX)), 
in the second derivative, d 2 e/d\ 2 , might be induced by K, 
inside that more trivial two-dimensional landscape corre- 
sponding to H. It must be kept in mind that H is often 
not very well suited to ensure a localisation of B's, a 
necessary condition for the exploration of an energy sur- 
face parametrized by 6's. A review of those landscapes 
obtained by constrained Hartree-Fock and/or Hartree- 
Bogoliubov, whether at zero or finite temperature, is in 
order, if only to analyze the role of fluctuations and also 
verify whether edges of variational spaces, and additional 
solutions, trapped by such edges, might exist. 
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